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Background: The next-to-next-to-next-to-leading order (N^LO) nuclear energy density functional extends the standard 
Skyrme functional with new terms depending on higher-order derivatives of densities, introduced to gain better pre- 
cision in the nuclear many-body calculations. A thorough study of the transformation properties of the functional with 
respect to different symmetries is required, as a step preliminary to the adjustment of the coupling constants. 

Purpose: Determine to which extent the presence of higher-order derivatives in the functional can be compatible with the 
continuity equation. In particular, to study the relations between the validity of the continuity equation and invariance 
of the functional under gauge transformations. 

Methods: Derive conditions for the validity of the continuity equation in the framework of time-dependent density functional 
theory. The conditions apply separately to the four spin-isospin channels of the one-body density matrix. 

Results: We obtained four sets of constraints on the coupling constants of the N'^LO energy density functional that guarantee 
the validity of the continuity equation in all spin-isospin channels. In particular, for the scalar-isoscalar channel, the 
constraints are the same as those resulting from imposing the standard U(l) local-gauge-invariance conditions. 

Conclusions: Validity of the continuity equation in the four spin-isospin channels is equivalent to the local-gauge invariance 
of the energy density functional. For vector and isovector channels, such validity requires the invariance of the functional 
under local rotations in the spin and isospin spaces. 

PACS numbers: 21.60.Jz, 11.30.-j 



I. INTRODUCTION 

In recent years, methods using energy density func- 
tionals (EDFs) [l| to describe nuclear properties are be- 
ing developed in three complementary directions. First, 
the ideas of effective theories 0, Q are einployed in deter- 
mining the EDFs from first principles 043 ■ These devel- 
opments are supplemented by a renewed interest [ sMlOl 
in the density-matrix expansion (DME) methods [liril2j|. 
which allow for treating exchange correlations in terms of 
(quasi)local functionals. Second, the coupling constants 
of the well-known EDFs undergo a thorough scrutiny, 
including an advanced work on the readjustment of pa- 
rameters [Tsl . IT^ and study of inter-parameter correla- 
tions . Finally, the standard functionals are extended 
by adding new terms so as to gain increased pre- 

cision of description and predictability, in_ quest for the 
spectroscopic-quality [IT 



21| EDFs. 



and universal 

In the present work we study properties of EDFs [16[ 
and pseudopotentials [l^ extended by adding terms that 
depend on higher-order derivatives up to sixth, next-to- 
next-to-next-to- leading order (N^LO). Such extensions 
lead to self-consistent mean-field Hamiltonians that are 
sixth-order differential operators [22^, that is, they de- 
pend on up to sixth power of the momentum operator. 
This makes them unusual objects, in the sense that stan- 



dard second-order one-body Hamiltonians contain only 
the Laplace operator in the kinetic-energy term and pos- 
sibly the angular-momentum operator in the spin-orbit 
term. The main question we address here is whether the 
presence of higher powers of momenta is compatible with 
the continuity equation (CE). 

The CE is a differential equation that describes a con- 
servative transport of some physical quantity (23j . In 
quantum mechanics, it relates the time variation of the 
probability density to the probability current . In our 
case, it appears when the N^LO EDFs or pseudopoten- 
tials are employed within a time-dependent theory. For 
the standard Skyrme (NLO) functional, the validity of 
the CE has been checked explicitly [25i]. Our goal here 
is to derive constraints on the coupling constant of the 
N'^LO EDF or parameters of the pseudopotential that 
would guarantee the validity of the CE. Apart from link- 
ing the CE to the local gauge symmetry [ll], we also 
analyze the CEs in vector and isovector channels and 
link them to the local non-abelian gauge symmetries. 



'Electronic address: | francesco.raimondi75@gmail.com] 



The paper is organized as follows. In Sec. Ill Al we 
present the standard quantal CE for a single particle and 
introduce the vector CE. Then, in Sec. IIIBI we discuss 
the CEs within the time-dependent density functional 
theory and in Sec. Ill CI we specify the case to the N'^LO 
quasilocal functional. The main body of results obtained 
for the CEs in the four spin-isospin channels is presented 
in Sec. |llT] and Appendices EUCl Finally in Sec. |IV] we 
formulate the conclusions of the present study. 
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II. CONTINUITY EQUATION IN THE EDF 
APPROACH 



A. Time evolution of a spin-^ particle 



We begin by recalling the well-known [2J| CE for a 
single particle. The time evolution of a non-relativistic 
spin-^ particle moving in a local potential is given by the 
Schrodinger equation, 

d 

ih—^{r<j,t) = - — Aij{ra,t) + Vo{r,t)i}{ra,t) 

+ J2 V^^{r,t)Y,{<yW^^W')i^ircT',t), (1) 



say, along the z axis, that is Vi{r,t) = Vi{r,t)ez- In 
this case, the time evolutions of the spin-up and spin- 
down components decouple from one another, that is, 
s(r,0) — s{r,0)ez implies s{r,t) = s{r,t)e2, and the 
spin-up and spin-down components individually obey the 
corresponding CEs. 

We also note here that for a nonlocal potential-energy 
term, 

(T>?A)(ra,t) = J d^r'^Vira,r'a',t)i;{r'a',t), (6) 

cr' 

the time evolution does not, in general, lead to a CE. 



where Vo{r,t) and Vi^(r,i) are scalar and vector real 
time-dependent potentials, respectively, and {a\<7^\a') 
are the standard Pauli matrices. By multiplying Eq. ^ 
with ip*(ra, t), summing up over cr, and taking the imag- 
inary part, we obtain the standard CE for the probability 
density p{r,t) in terms of the current j{r,t), 



d_ 
di 



m 



-V-j{r,t), 



where 



p{r,t) - ^|V(ra,i)P, 



(2) 

(3a) 
(3b) 



We see that the hermiticity of the local potential guar- 
antees that the potential energy does not contribute to 
the CE of Eq. ©. 

Similarly, by multiplying Eq. ([T]) with 
'ijj*{ra" ,t){a"\a^\a), summing up over a" and cr, 
and taking the imaginary part, we obtain the CE for the 
spin density s^{r, t) in terms of the spin current Ju{f, t), 

^Ur,t) = --V • Mr,t) + l(Viir,t)xs{r,t)) ,(4) 
ot m n\ J V 



where 



(5a) 



J^(r, t) = Y^ Im( V/(ra', t)(a'|a,|a) VV'lra, t) ) . (5b) 



We see that the spin CE does depend on the vector po- 
tential, and the second term in Eq. Q is responsible, 
e.g., for the spin precession in magnetic field. 

It is interesting to note that when potential Vi(r,i) is 
parallel to the spin density s(r, t) (non-linear Schrodinger 
equation), all components of the spin density fulfill the 
CEs. In fact, this is exactly the case for the TDHF equa- 
tion induced by a zero-range two-body interaction, see 
below. Another interesting case corresponds to the vec- 
tor potential aligned along a fixed direction in space. 



B. Time-dependent density functional theory 

In the framework of the time-dependent Hartree-Fock 
(TDHF) approximation or time-dependent density func- 
tional theory (TDDFT), the so-called memory effects are 
often neglected and it is assumed that the potential at 
time t is just the static potential evaluated at the in- 
stantaneous density |27|. For these two time-dependent 
approaches, the starting point is the equation of motion 
for the one-body density matrix pap [H, [1^ , 



(7) 



where the mean-field Hamiltonian hap is defined as the 
derivative of the total energy E^p] with respect to the 
density matrix. 



dE{p\ 

dpjia 



(8) 



In the present study we are concerned with the Kohn- 
Sham approach [29j , whereby the total energy is the sum 
of the kinetic and potential-energy terms. 



where 



E{p) = Ek{p'\ + E^{p\, 



Ek{p} = ^Jd'rrS{r,t) 



(9) 



(10) 



and ro°(r,t) = (E^tt'^ ' 'V' p{raT,r'aT,t)j\r=r' is the 
scalar-isoscalar kinetic density, see, e.g., Ref. [s^l for def- 
initions. The nonlocal density, can be defined in terms 
of either the fixed-basis orbitals, ipa{fO'T), 

p{raT,r'a'T',t)=J2Mr''^)Ppo.{t)rjr'a'T'), (11) 



or instantaneous Kohn-Sham orbitals, (j)i{raT,t), 

A 

p{raT,r'a'T',t) ^^(l),{raT,t)(t)*{r'a'T',t). 

i=l 



(12) 
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The mean-field Hamiltonian is the sum of kinetic and 
potential-energy terms, h^is — T^ib + T^p, where 



and 



dEp{p} 

dp 13 a 



(13) 



(14) 



Let us now assume that the potential energy is invari- 
ant with respect to a unitary transformation of the den- 
sity matrix (26l.[28j. U = exp{i?]G), that is, for all param- 
eters T] we have, 



Ep{p} = Ej,{UpU+}, 



(15) 



where Gais is the hermitian matrix of a one-body sym- 
metry generator. Then, the first-order expansion in rj, 



Ep{UpU+}^E,{p} + rjJ2 

l3a 



dEp{p} d{UpU+)p^ 



dpi3a 



drj 



7-1 = 

... (1^) 

gives a condition for the energy to be invariant with re- 
spect to this unitary transformation, that is 



Trr[G,p] = TrG[r,p] = 0, 



(17) 



which allows us to derive the equation of motion for the 
average value of (G) = TvGp. Indeed, from the TDDFT 
equation (O we then have: 

ih^ (G) = ihTiG^p = TrG[/i, p] = TrGlT, p] , (18) 
dt at 

that is, the time evolution of (G) is governed solely by 
the kinetic term of the mean-field Hamiltonian. 



1. Continuity equation for the scalar-isoscalar density 

The CE now results from specifying r^G to the local 
gauge transformation (2^, |3l| that is defined as 



ip'^{raT) = {Utl;a){r<7T) 



— pi-fir) 



(19) 



Then, Eq. ^ gives: 

p'irar, r'a'r', t) = e^'^^''^-'^^'''^) p{raT, r'a'r', t). (20) 

Matrix elements of the local-gauge angle 7(r) are given 
by local integrals. 



7a/3 



j d^r ip'^{raT)^{r)iljfj{raT) 



(21) 



therefore, from Eq. (ITT|) again, the average value of the 
gauge angle, (7) — Trjp, depends on the scalar-isoscalar 
local density p°{r, t) = J2aT Pi'''^"'', far, t), that is. 



(7) - / d3r7(r)p°(r,0. 



(22) 



Now, the assumed local-gauge invariance of the poten- 
tial energy implies the equation of motion for the average 
value (7), which from Eq. (fT8| reads 



dt 



(7) 



h 

m 



(23) 



where the standard scalar-isoscalar current is defined as 
M foir,t) = Ep^M^ - V')/o(rar,r'aT,i)],^,,. 

We note here [26ll31|. that the gauge invariance that 
corresponds to a specific dependence of the gauge angle 
on position, -f{r) = Pq ■ r, represents the Galilean in- 
variance of the potential energy for the system boosted 
to momentum Pq. Then, equation of motion p3p simply 
represents the classical equation for the center-of-mass 
velocity, 

d (r) _ d (P) _ {-^hV) 

di^ = d^^^^ ^^^^ 

In the general case, that is, when the potential energy 
is gauge-invariant and the gauge angle 7(r) is an arbi- 
trary function of r, Eq. ([23| gives the CE that reads 



i.S(M)^-iv..-S(M)^ 



(25) 



Thus for a gauge-invariant potential energy density, the 
TDHF or TDDFT equation of motion implies the CE, 
that is, the gauge invariance is a sufficient condition for 
the validity of the CE. By proceeding in the opposite di- 
rection, we can prove that it is also a necessary condition. 
Indeed, the CE of Eq. implies the first-order condi- 
tion (jl7p . and then the full gauge invariance stems from 
the fact that the gauge transformations form local U(l) 
groups. 



Continuity equation for densities tn sptn-isospm 
channels 



We can now repeat derivations presented in Eqs. ([TJ 
by considering the spin-isospin local-gauge groups, 
and derive CEs in other spin-isospin channels. To this 
end, we first express the nuclear one-body density matrix 
(HU^dl) as a linear combination of nonlocal spin-isospin 
densities p^(r,r') [9], 



p{raT, r'a'r'] 



v=0,l 
t=0,l 



V3 



v+t 



,,p*(r,r')]1 , (26) 



where the sums run over the spin (v = 0,1) and isospin 
(t = 0, 1) indices denoted by subscripts and superscripts, 
respectively, coupled to total scalar and isoscalar. Here 
and below we use the coupling of spherical tensors both 
for angular momentum and isospin tensors; therefore, in 

Eq. jMl) the factor of (v^)" was included so as to can- 
cel the corresponding values of the Clebsch-Gordan co- 
efficients, and to maintain the standard normalization 
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of the spin-isospin densities. The spin-isospin densities 
can be conversely expressed as the foUowing traces of the 
density matrix, 

Pl{r,r') ^ ^'^r'r^rPirar^r'aW'). (27) 



The CEs for densities in the scalar-isoscalar {v = 0, 
t = 0), scalar-isovector [v = 0, t = 1), vector-isoscalar 
{v — 1, t = 0), and vector-isovector = 1, t = 1) 
channels, 



Ap*(r.) = -Av.J*(r), 

at m 



(28) 



where J*(r) = j-^ (V ~ V) pl{r,r')\r'=r and /9*(r) ^ 
Py{r,r), are now equivalent to the local gauge invari- 
ances, respectively, with respect to the four local spin- 
isospin groups: 



(29) 



Of course, the standard CE derived in Sec. IIIB II cor- 
responds to 7(r) = 7o(r). Note that the four gauge 
groups are different: J7o(r) gives the standard abelian 
gauge group U(l), C/i (r) and U^{r) form the non-abelian 
gauge groups SU(2), whereas Ul{r) corresponds to the 
non-abehan gauge group SU(2)xSU(2). 



C. The N^LO quasilocal functional 

We are now in a position to discuss the CE for the 
N'^LO quasilocal functional introduced by Carlsson et 
al. ^16i] . By imposing on the functional the gauge- 
invariance conditions, we can then confirm and explicitly 
rederive the results of Sec. IIIBI The explicit derivation 
will also allow us to discuss the CEs for densities in other 
spin-isospin channels analyzed in Sec. Ill B 21 

Below we consider the EDF given in terms of a local 
integral of the energy density 'He{i'), 



E{p} = J d^r HEir), 



(30) 



which is represented as a sum of the kinetic and potential 
energies conforming to Eq. 



(31) 



t=oa 



To lighten the notation and avoid confusion with the 
isospin index i = 0, 1, in this section we do not explicitly 
show the time argument of densities, which within the 
TDDFT ah depend on time. 

The quasilocal N'^LO EDF was constructed JJi] by 
building the t — Q and t — 1 potential-energy densities 
H*{r) from isoscalar and isovector densities, respectively, 
and their derivatives up to sixth order. For clarity, we 



give here a brief summary of definitions and notations 
used in this construction. 

The local higher-order primary densities are defined 
by the coupling of relative-momentum tensors KnL [3] 
with nonlocal densities (j27l) to total angular momentum 
J, that is. 



(32) 



Then, a general term of the N'^LO functional can be writ- 
ten, in the language of the spherical tensors, as 



rpn'L'v' J' .t, •._ 
^mI,nLv,J \' )- 



[pn'L'.'j'{r) [DraipUvAr)]^], (33) 



where the local secondary densities, [^m/P^L„j('")] j/, are 
obtained by acting with derivatives D^j on primary den- 
sities and coupling them to total angular momentum J'. 
Each term (I33p is multiplied by the corresponding cou- 



pling constant C^j^^j^f that is denoted by the same set 
of indices as those in the term itself. 

We note here that the definition of the isovector terms 
depends on whether one uses Cartesian or spherical rep- 
resentation of tensors in isospace. On the one hand, the 
use of the standard Cartesian representation, see, e.g., 
Rcfs. [16, 30], implies that the isovector terms depend on 
products of differences of neutron and proton densities. 
On the other hand, the use of the spherical representa- 
tion, which was assumed in Ref. [l^ and is also used in 
the present study, involves the coupling of two isovectors 
to a scalar, whereby there appears a Clebsch-Gordan co- 
efficient of (a/S) . Therefore, for the isospace spherical 
representation, the isovector coupling constants are by 
the factor of (a/S) larger than those for the Cartesian 
representation. 

In the remaining part of this sectioUj we employ 

the compact notation introduced in Ref. [22| . whereby 
the grouped indices, such as the Greek indices a = 
{riaLaVaJa} and the Roman indices a = {mala}^ denote 
all the quantum numbers of the local densities Pair) and 
derivative operators Z?a, respectively. In this notation, 
the N^LO potential-energy density of Eq. reads 



(34) 



aal3 



Our following discussion of the CE is mainly focused on 
the one-body potential-energy term, defined in Eq. 
as the variation of the potential energy with respect to 
the density matrix. For the N'^LO functional, this term 
was derived in Ref. [22| , where it was shown that in space 
coordinates it has the form of a one-body pseudopoten- 
tial. 



Kir) 



(35) 



An equivalent form of the one-body pseudopotential, 
which can be obtained by recoupling spherical tensors 
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within a scalar, and which separates out the spin Pauh 
matrices, reads 



(36) 



In turn, potentials C/^(r) were derived as linear com- 
binations of the secondary densities. 



aaj3;d5 



(37) 



7^ 

where Xa'a% numerical coefficients. We call the one- 
body operator F^^, (r) pseudopotential, because it is de- 
fined in terms of potentials U^{r) and differential opera- 
tors Dn^L^ acting on single-particle wave functions. We 
note here, that in Eqs. ([55)) and ([55)) potentials always 
appear to the left of all derivatives; nonetheless, the one- 
body pseudopotential is a hermitian operator, which is 
guaranteed by specific conditions obeyed by potentials 
?7^(r), which were derived in Ref. (23 |. 

For the one-body pseudopotential ([35)) . the 
Schrodinger equation that gives the time evolution 
of single-particle Kohn-Sham wave functions in space 
coordinates reads, 

d 

ih—(j)i{raT,t) = -—A4)i{r(TT,t) 

+ Y.Kr:ir)Mr<j'r',t). (38) 

cr'r' 

By multiplying the Schrodinger equation with the 
complex-conjugated wave function, (f>*{r'a'T' ,t) and 
summing over the single-particle index i we obtain the 
time-evolution equation of the density matrix ()12p , that 
is, 

d 

ih—p{raT,r'a'T't) = 

(A - A') p{raT, r'a'r', t) 

2m 

+ J2 {KrAr)p{ra"T",r'a'T',t) 

a" r" 

- vi.%*(r')p{raT, r'a"T\ t)^ . (39) 

Before we proceed, we must first consider the complex- 
conjugated pseudopotential r^,^„*(r'). To this end, we 
use the property of the Biedenharn-Rose phase conven- 
tion employed in Refs. [16„i22], by which all scalars are 
always real. Note that for the spherical representation of 
Pauli matrices, the Biedenharn-Rose phase convention 
implies the transposition of spin indices, that is. 



(<;')* = (-ir 



(40) 



where fi = for v = and fi = —1, 0, 1 for v = 1 denote 
tensor components of scalar and vector Pauli matrices, 
respectively. 

Finally, in Eqs. ([55)) and (1571) . the complex 
conjugation only affects coefficients xlfaf-f = 



f^;?;;*(r') = fXT'(r'), 



(41) 



for 



f:-'(r') = E 
and 



[c/;*(r')^;^i^ 



(42) 



aa/3;d(5 

(43) 

It means that in all further derivations we must use the 
second set of potentials U^{r') with signs of terms mod- 
ified according to the phase (— l)"7+"'<'+'"ti. It is now 
obvious that the CEs will hold independently of the spin- 
isospin coordinates if, and only if, the pseudopotentials 
fulfill the condition 



E (r?,<(r)p(ra"r",r'a'r') 



FXT'(r')p(rar,rV"r") 



0. (44) 



We are now in a position to separate the four spin- 
isospin channels in Ec^. ([39)) . We do so by multiplying 
both sides of the equation with ct^ '^'''t't ^'^'^ summing 
over aT^a'r'. From Eq. ([27| it is then obvious that, 
in close analogy to Sec. Ill A[ after setting r' = r, we 
obtain the CEs ([28)) in the four spin-isospin channels, 
provided terms coming from one-body pseudopotentials 
do not contribute, as in Eq. ([44)) . When evaluating this 
condition for the four spin-isospin channels, we use the 
expression for the trace of three Paul matrices in spheri- 
cal representation, which reads [s^j, 



Tr 



+v"){-ir 

(45) 

where we introduced A{v+v' +v") as a shorthand symbol 
for numerical coefficients coming from the computation 
of the trace. In the calculation, we only need values of 
yl(0) = 2, A{2) 2^/3, yl(3) = 2\/2i. After a trivial but 
lengthy calculation, we obtain the final result: 
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y A{v + v' + v^)A{t + t' + t")^ -r, ^{2J+1)\ ' 

7,t"J 



J 



(r) A.,L,P*',(r,r') 







J. 


V 



+Vy-V ^_^Y' +t" ~t 
I 



U'* (r') D'j:,{r,r') 



0. (46) 



For our practical implementation of the CE condi- 
tion (|46p , we proceed by transforming the two differential 
operators, Dn^^^ and D'^^j^^, which act on two differ- 
ent variables r and r', respectively, with the recoupling 
methods developed in Ref. [22j, and we obtain, 

Dn^Ly = E C2^7^,(*)"2-"'[A.L'i^„L]i^,(47) 



^Ln' L' 



nLn' L' 



On the right-hand sides, operators K^l are the higher- 
order spherical tensor derivatives 18| built of the relative 
momenta, k = (V — V')/2i, and operators Dn'L' act 
on variable r after one sets r' = r. The 91 numerical 
coefficients K^J^j^,, which are needed up to N'^LO, have 
been derived in Ref. [22*1. By using Eqs. (gZl) and (l48t . 
one can express the last line of the Eq. (|46p as a linear 
combination of products of pairs of secondary densities 
coupled in the spin and isospin spaces to ranks v and 
t, respectively. This final form, which for brevity is not 
shown here explicitly, is used in obtaining the results of 

secnni 



ones through independent ones. Obviously, for a given 
set of linear conditions, this can be done in very many 
different ways; below we present in each case only one 
choice thereof. We also use the name of a free coupling 
constant to denote either the unrestricted or independent 
one. 

The structure of this section is as follows. First, in 
Table (HI we present for the four spin-isospin channels an 
overview of results by showing the number of unrestricted 
(U), vanishing (V), independent (I), and dependent (D) 
coupling constants of the EDF. The sum of the unre- 
stricted and independent coupling constants gives the 
number of the free ones. Second, in Sec. IIII Al we dis- 
cuss the simplest case of the zero-order terms, where one- 
body pseudopotentials reduce to simple potential func- 
tions, and we can link our results to those presented in 
the introductory Sec. Ill Al Next, in Sec. IIII B I we briefly 
describe the results obtained for the standard CE in the 
scalar-isoscalar channel. For the CEs in the three other 
spin-isospin channel, we present our results in Secs. lIIl"Cl - 
IIIIEI where, for clarity, only the second-order terms and 
general rules are discussed, while the results for fourth 
and sixth orders are collected in Appendices [SHCl 



III. CONTINUITY EQUATIONS IN THE FOUR 
SPIN-ISOSPIN CHANNELS 

Condition sets constraints on the coupling con- 
stants Cf'* of the EDF. In our study, these explicit 
constraints were obtained, with the aid of the symbolic 
programming, as solutions of a linear system of equa- 
tions, where each equation is found by considering the 
coefficients standing at a given product of pairs of sec- 
ondary densities. Indices correspond to the choice 
of the channel under examination: (0,0) for the CE in the 
scalar-isoscalar channel, (0,1) for the CE in the scalar- 
isovector channel, (1,0) for the CE in the vector-isoscalar 
channel, and (1,1) for the CE in the vector-isovector 
channel. 

As a result of this analysis, for each spin-isospin chan- 
nel we can classify the coupling constants in four cate- 
gories defined in Ref. [Tcj . namely, unrestricted, vanish- 
ing, independent, and dependent. The unrestricted cou- 
pling constants are not affected by condition p6)) and 
vanishing ones are forced by this condition to be equal 
to zero. The remaining coupling constants obey sets of 
linear conditions, whereby one can express the dependent 



A. Constraints for zero-order terms 

For zero-order terms, condition (|46|) gives no con- 
straints on the coupling constants, apart from those in 
the vector-isovector channel, where we have, 

^0000,1 _ 1 p0011,l 

'-'00,0000 ~ ^^2^-^00,0011' l^'iaa; 

^0011,0 _ 1 ^0011,1 /.qK^ 
'-'oo.ooii — ^^2^-^00,0011' yi-auj 

whereas the coupling constant Cqq po'oo unrestricted. 

It is interesting to discuss these results in connection 
with the derivation of the CE for a spin- ^ particle moving 
in a local potential, which we gave in Sec. Ill Al There, 
we pointed out that the CE in the vector channel is 
valid when the vector potential is parallel to the spin 
density. Exactly this situation occurs for the zero-order 
EDF, where in each spin-isospin channel the potential 
functions are simply proportional to densities. 

In fact, the simple algebraic rule of the vector prod- 
uct in Eq. Q is equivalent to the coupling of pairs of 
identical commuting rank 1 tensors to rank 1, which is 
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TABLE I: Number of unrestricted (U), vanishing (V), independent (I), and dependent (D) coupling constants of different orders 
in tlie EDF up to N'^LO, sliown for the four spin-isospin channels. 



Order 


Total 




v = 0, t 


= 






v = 0, t 


= 1 






V =1, t 


= 






V ^l, t 


= 1 




U 


V 


I 


D 


U 


V 


I 


D 


U 


V 


I 


D 


U 


V 


I 


D 





4 


4 











4 











4 











1 





1 


2 


2 


24 


6 





8 


10 


3 


7 


6 


8 


2 


10 


4 


8 


1 


11 


1 


11 


4 


90 


6 


54 


6 


24 


3 


57 


6 


24 


2 


64 


4 


20 


1 


65 


1 


23 


6 


258 


6 


200 


6 


46 


3 


203 


6 


46 


2 


216 


4 


36 


1 


217 


1 


39 


N^LO 


376 


22 


254 


20 


80 


13 


267 


18 


78 


10 


290 


12 


64 


4 


293 


4 


75 



identically null. This is just the case for the coupling to 
rank v = 1 (t = 1) in the spin (isospin) space for the 
vector-isoscalar (scalar-isovector) channel of the CE, in 
which the identically null tensors formed by the pairs of 
densities at zero order leave the corresponding coupling 
constants unrestricted. However, when the coupling to a 
rank-1 tensor is simultaneously performed in both spin 
and isospin space (v = 1 and t = 1), as is the case in 
the vector-isovector channel of the CE, the two negative 
signs in the commutation of the pair of densities give an 
overall positive sign, and the above selection rule does 
not apply. This explains why, at zero order of the vector- 
isovector channel, condition (|46p does induce constraints 
on coupling constants - those given in Eqs. (|49ap and 
(l49b)) . 



the functional. At second order we obtain, 



0000, t 
00,2000 



a 



00,1110 



C, 



00,1111 



a 



1112,* 
00,1112 



^0000, t 
'-'ll.llll 



c: 



0000.1 
20,0000 



3* ^1101,1- 

^^00,1101 



3^3 
3* 



r' 



(50a) 



00,2011 2 Y 2^00,2211 

'(50b) 



^0011,1-t J /T^0011,l-i (rr. . 

2 '-^00,2011 +"^VOOqq 2211 , [OVC) 



— — / ^/^0011.1-t _ _3_ 
3 V 3 °"^20ii 

^1101,* _ 
•-^11,0011 ~ 

^0011.1 _ ^0011,1 _ p, 
•-^20,0011 ~ "-^22,0011 ^ 



^0011, 1-t fr:nA\ 
^nn 2211 i'^'-J^) 



(50e) 
(50f) 



B. Constraints for the scalar-isoscalar channel 

For the standard CE in the scalar-isoscalar channel, 
Eq. (pS)) . the constraints that we found from Eq. are 
exactly the same as those defining the gauge-invariant 
functional up to N^LO This constitutes an explicit 
verification of the general result presented in Sec. Ill B II 
As in Ref. [11] , we obtained two sets of independent con- 
straints, one for the isoscalar coupling constants Cf '° and 
one for the isovector coupling constants Cf'^. Within 
each isospin channel, the linear conditions for the scalar 

coupling constants, ii'o'i! ' ^'"^ disconnected from 

the linear conditions for the vector coupling constants, 

c: 



npLp IJff ,t 



Constraints in Eqs. (|50ap - (|50dp connect the isoscalar 
and isovector coupling constants. The numerical coeffi- 
cients of the corresponding linear combinations are the 
same as those for the scalar-isoscalar CE, see Eqs. (Cl)- 
(C4) of Ref. [ll], apart from factors of -s/S explained 
before Eq. (p4l) . However, conditions for the scalar- 
isoscalar CE keep the coupling constants in the two 
isospin channels disconnected. Moreover, in the scalar- 
isovector channel the spin-orbit coupling constants must 
vanish, Eq. (|50ep . along with the isovector surface cou- 
pling constants, Eq. (j50f|) . On the other hand, the cor- 
responding isoscalar surface coupling constants, C'2o,odooi 
C200011, and ^22 00115 are left unrestricted. 

For the fourth and sixth orders, analogous constraints 
are presented in Appendix [Xj 



Constraints for the scalar-isovector channel 



D. Constraints for the vector-isoscalar channel 



Validity of the CE for the scalar-isovector density, 
Eq. (pS)) for D = and t — 1, imposes through Eq. pS)) 
certain specific constraints on the coupling constants of 



Validity of the CE for the vector-isoscalar density, 
Eq. for t; = 1 and t = 0, imposes through Eq. pS)) 
at second order the following constraints on the coupling 
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constants of the functional, 



^1101, t 
'-"00,1101 



1110, t 
00,1110 



-—C 



,0011.* 



,0000.* 



0000,* 
00,2000' 



'J ^0000.* 
2 '-"00,2000' 



0011,* 
20,0011 
^0000,* 

'-^11,1111 

0011,* 
00,2211 



c: 



•-"22,0011 

^1101,* 

•-"11,0011 



0, 
0, 



(51a) 

(51b) 
(51c) 

(51d) 

(51e) 
(51f) 
(51g) 



whereas the two coupHng constants ^290000 l^ft unre- 
stricted. We note here that the constraints now connect 
scalar and vector coupling constants. Altogether, at sec- 
ond order, for the vector-isoscalar channel of the CE we 
have 6 free and 8 dependent coupling constants. Apart 
from that, 10 second-order coupling constants must van- 
ish, which includes the surface ones in Eq. (j51el) . spin- 
orbit ones of the Eq. (j51fl , and tensor ones in Eq. ( |51g[ ) . 

For the fourth and sixth orders, analogous constraints 
are presented in Appendix |B] 



E. Constraints for the vector-isovector channel 

Validity of the CE for the vector-isovector density, 
Eq. for w = 1 and t = 1, imposes through Eq. pS)) 
constraints on the coupling constants that relate them in 
both spin and isospin spaces. At all the orders, we can 
express all dependent coupling constants through only 
one vector coupling constant, which can be chosen either 
in the set of the vector-isoscalar or vector-isovector ones. 
At second order, these constraints read. 



^0000.* 

•-"00,2000 




(52a) 


^1101.* 
•-"00,1101 




(52b) 


^1110,1 
•-"00,1110 


— vooqoiiio, 


(52c) 


^1111,* 

•-^00,1111 




(52d) 


^1112,* 
•-"00,1112 




(52e) 


^0011.* 
•-'00,2011 


— 0(\ 0) L^00,1110' 


(52f) 


^0011,* 
•-"20,0011 


^0011,* p. 
— "-"22,0011 ~ 'J' 


(52g) 


^0000.1 
•-"20,0000 


= 0, 


(52h) 


^0011,* 
•-"00,2211 


= 0, 


(52i) 


^0000,* 
•-"11,1111 


^1101,* 
~ '-'11,0011 ~ '-'' 


(52j) 



second-order coupling constants, which includes the sur- 
face ones in Eqs. ( |52gP ~ (j52h|) . the tensor ones in Eq. (|52i|) 
and spin-orbit ones in Eq. ( |52j[ ), are forced to be equal 
to zero. 

For the fourth and sixth orders, analogous constraints 
are presented in Appendix O 



IV. CONCLUSIONS 

In the present work, we have derived sets of constraints 
on the coupling constants of the N^LO energy density 
functional that guarantee the validity of the continuity 
equation in the four spin-isospin channels. In the scalar- 
isoscalar channel, these constraints are identical to those 
induced by the standard local gauge invariance condi- 
tions. We extended this connection to vector and isovec- 
tor channels, where the validity of the continuity equa- 
tions is equivalent to the local gauge invariance with re- 
spect to spin and isospin rotations, respectively. 

We note here that in our analysis we implicitly as- 
sumed that all densities that build the N'^LO energy 
density functional are nonzero. Of course, there can be 
many specific situations when some densities vanish, and 
thus the coupling constants related to them become un- 
restricted. This occurs, for instance, when densities are 
restricted by some symmetry conditions. Obviously, the 
methods developed in the present work can then be ap- 
plied to derive new (weaker) sets of constraints that cor- 
respond to each one particular case. For example, when 
the proton- neutron symmetry is conserved, see, e.g., 
Refs. [13, HI], the ±1 components of all isovector densities 
{v = 1) vanish, and the energy density is invariant with 
respect to a one-dimensional U(l) gauge rotation in the 
isospin space, that is, with respect to U = exp(i7ioT^°). 
Then, the continuity equations for neutrons and protons 
decouple from one another and become independently 
valid, provided the isoscalar and isovector coupling con- 
stants independently obey the standard gauge-invariance 
conditions. 
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Appendix A: Constraints for the scalar-isovector 
channel (fourth and sixth orders) 



and only coupling constant C'. 



0000,0 
20,0000 



is unrestricted. Al- 



together, at second order for the vector-isovector channel 
of the CE we have 2 free coupling constants and 11 cou- 
pling constants that are dependent. All the remaining 11 



At fourth order we found that the isoscalar and isovec- 
tor coupling constants are connected or not depending 
on the parity of quantum numbers n. Moreover, simi- 
larly as for the scalar-isoscalar channel of the CE, the 
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scalar and vector coupling constants are kept apart. The 
constraints among the scalar coupling constants read, 



r' 



0000, t 
00,4000 

2000, t 
00,2000 



0^2202, t 
■J '-'00, 2202 

1 ^/F/^2202,t 

2 ^ ^"-^00,2202' 



^1101, t _ _r, /3 „t „2202,l-t 
'-^00,3101 ~ g"^ '-'00,2202 ; 



(Al) 
(A2) 
(A3) 



and those among the vector coupling constants read, 



C, 



1110. t 
00,3110 



_^otr'2212,l-t _ _^qt^0011,l-t j 
' ^/g '-'00,2212 '-^00,4211 A^^J 



^oo,3i'ii — ^■^(3*)y'^c', 



2212, l-t 
00,2212 ' 



(A5) 



00,3112 

1112.* 
00,3312 

'0011, t 
00,4011 

2011. t 
00,2011 

2011. t 



r,/ot^i/-<2212,l-t _ 14 3* „0011.1-t ( 
^l'^ J'-'00,2212 g _y2'-^00,4211 ' \^^) 



2 /'''^0011,1-t 
'3l_ty 5^00,4211 , 

nr 7^0011, t 

-J / -J ^2212,t . ' '-'00,4211 



2 V 5 

1 

2 



00.2212 



-t- 



4V5 



(A7) 
(A8) 



^2211, t _ I " ^. 
L'nn 9211 ^ \/ g '-'00,2212 



'-'00,2213 ~ \/ g'-'00,2212 "T 2 y "5" 



15CoO,2212 + Y2^^0M211' (^^) 
rC'nn ,19111 (-^10) 

(All) 
(A12) 



•J -^2212, t , "-'00,4211 



3V5 ' 



7 „2212,t , 1/21^0011,* 



We also found that the fourth-order surface isovector 
coupling constants must vanish. 



0000,1 

40,0000 



^0011,1 _ puuii,i _ p, /'A^Q^ 
'-^40,0011 ^ '-^42,0011 ~ '^' l^-'-'Jj 



^0011,1 



whereas the corresponding isoscalar coupling constants 
are unrestricted. Apart from the coupling constants dis- 
cussed above, all the remaining 54 fourth-order coupling 
constants are forced to be equal to zero. 

In the same way, at sixth order we found the following 
constraints for the scalar. 



^0000,* 
'-^00,6000 ~ 

^2000, t 
'-'00,4000 


0^3303, l-t 
2* "^'-^00,3303 

aV7 ' 

0*3/^^3303,1-4 
4 '^00,3303 ' 


(A14) 
(A15) 


^2202, t 
'-^00,4202 — 


otq / 5 ^3303, l-t 
V 7 '"''3303 ' 


(A16) 


^1101,* 
'-'00,5101 — 


9 /3„3303,t 
2 V 7 '"''3303' 


(A17) 


^3101,* 
'-'00,3101 — 


9 /9Tr'3303,t 
^q"^ ^J^'-' 00, 3303' 


(A18) 



and vector coupling constants. 



00,5110 

1111, t 

00,5111 

1112, t 
00,5112 

,1112, t 
00,5312 

3110, t 
00,3110 

3111, t 
00,3111 

3112, t 
00,3112 

3112, t 
00,3312 

3312, t 
00,3312 

3313,4 
00,3313 

3314, t 
00,3314 

,0011, t 
00,6011 

2011, t 
00,4011 

2011, t 
00,4211 

2211, t 
00,4011 

2211, t 
00,4211 

2213, t 
00,4213 

2213, t 
00,4413 



^2212. l-t 
,4 '-'00,4212 



3*3 



3^0011, i-t 
5^" 



-'00,6211 



~,2212.1-t 

■'00.4212 ' 



t 1 ^2212, l-t 



^ 2^' 



qt 6 /q^0011,l-t 

5 00,6211 



31-* V 5 
-^2212, i-t 

gt ''-'00,4212 



r0011,l-t 

00,6211 ' 



(A19) 

(A20) 

(A21) 
(A22) 



>21 



IOV5 



'£^0011,1-* A2q^ 
[^'-^00,6211 lA^'Jj 
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C, 



00,4212 



(A24) 



Qt ^2212, l-t 
10 '"^^^^ 



,42^ 
25' 



3* — \/3Cqq g2ii (A25) 



12 7 



^00,6211 



(A26) 



_qt ^2212.1- 

g '-'00,4212 

_3^ /7„2212,l-t 
9 V 5 00,4212 ' 



3*^^C^^iY*' (A27) 



2212, i-t 



8C 



i0011,l-t 
00.6211 



*V5 



1 

4 

7 

4 

6Ci 
21 



3C 



Vl5 

'0011, t 



-"00,4212 



2^5 ' 
0011, t 



00.6211 



21c; 



5 

0011, t 

00,6211' 



00,4212 



2V5 



' r' 



00,6211' 



3 „2212,t 
g '-^00,4212 



12c 



00,6211 



0, 



2212, t 
00,4212 



18 



3 



C 



flOll.f 



4C 



0011, t 

00,6211 



V5 



(A28) 

(A29) 

(A30) 

(A31) 
(A32) 
(A33) 

(A34) 

(A35) 

(A36) 



We also found that the sixth-order surface isovector 
coupling constants must vanish. 



r' 



0000,1 

60,0000 



r 



60,0011 



r 



62,0011 



0, 



(A37) 



whereas the corresponding isoscalar coupling constants 
are unrestricted. Apart from the coupling constants dis- 
cussed above, all the remaining 200 sixth-order coupling 
constants are forced to be equal to zero. 

We have seen that for the scalar-isovector channel, the 
coupling constants are diagonal or nondiagonal in the 
isospin quantum number t. We can understand this point 
considering the fact that in order to separate the scalar- 
isovector channel of the CE and obtain condition for 
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t = 1, we have to multiply Eq. (|39| by the isospin op- 
erator Tt't- Then, the isospin index t" tells us in which 
half of the isospin space the coupling constants is. Non- 
diagonal constraints mean, in fact, that the same pair 
of secondary densities in the final form of condition p6)) 
can be produced by two terms of the functional that are 
isoscalar and isovector. This is possible, because the cou- 
pling to rank t = 1 allows for pairs of densities nondiag- 
onal in the isospin space, and this, in turn, boils down 
to constraints for coupling constants nondiagonal in the 
isospin space. 



Appendix B: Constraints for the vector-isoscalar 
channel (fourth and sixth orders) 

At fourth order, we again found two identical sets of 
linear combinations of the isoscalar and isovector cou- 
pling constants, in which the scalar and vector coupling 
constants are in same case connected to one another, that 
is. 



and we have, 



0000, t 
00,4000 



3C, 



00.2202 



2\/5 



2000,* 
00.2000 



-'- /F/-<2202,t 



6C, 



00,2202' 
2212, t 



00,3101 



a 



1110,* 
00,3110 



= -2 



C, 



6C, 



V5 ' 

3 „2202,t 
^•-^00,2202' 

2202,* 



00,3111 



C, 



00,3112 



V5 ' 

r, A7/-'2202,* 
-ZV OL/QQ 2202; 



0011,* 
00,4011 



00,2011 



3^2212,* 
g'-'00,2212; 



2212,* 
00,2212! 



2211,* 
00,2211 



2213,* 
00,2213 



2011,* 



3^2212,* 



7 „2212,t 
_ *^nn 991 9 



^1112,* _ ^0011,* 
•-^00,3312 — '-^00,4211 



= 0, 



(Bl) 
(B2) 
(B3) 
(B4) 

(B5) 
(B6) 
(B7) 
(B8) 
(B9) 

(BIO) 
(Bll) 



'0000,* 
40,0000 



left unrestricted. 



with the two coupling constants C, 
Apart from these 6 free and 20 dependent coupling con- 
stants, the vector-isoscalar channel of the CE requires 
that all the remaining fourth-order coupling constants are 
forced to be equal to zero. In particular in the Eq. (|B11[) 
we showed the vanishing coupling constants, belonging to 
the set of ones with indices m = and / = 0, which were 
found to be non- vanishing in the scalar-isoscalar channel. 
At sixth order the pattern of the results is the same. 



r" 



3303,* 
00,3303 



r' 



oooo,* 

00,6000 



■2000,* 
00,4000 



•2202,* 
00,4202 



1101,* 
00,5101 



3101,* 
00,3101 



C. 



1110,* 
00,5110 



1111,* 

00,5111 



^3110,* 
•-"00,3110 



3111, * 
00,3111 

3112, * 
00,3112 



^3312,* 
•-"00,3312 



'^1 



3313,* 
00,3313 



r" 



■3314,* 
00,3314 



0011,* 
00,6011 



00,4011 



^2213,* 
•-"00,4213 

^2211,* 
•-"00,4011 
^3112,* 
•-"00,3312 



15 

1112.* 
00.5112 



00,4212' 



2V15 

1112,* 
00,5112 



7C, 



2V15 



2C, 



1112,* 
00,5112 



3C, 



V3 

2212,* 



00,4212 

"271"' 

2212,* 



21G 



10\/5 



1112,* 
00,5112 



3^1112,* 
g'-'OO, 51121 



7C, 



1112,* 
00,5112 



5^5 



/3„1112,* 

sV 5 "^O'^"'' 



7 



,1112,* 



^•-^00,5112' 

2 ^1112,* 
g^-"00,5112J 

2 /7„iii2,* 



"00,51121 



2C, 



1112,* 
00,5112 



3^/5 ' 

4Y ^•-^00,4212: 
Z /^/-•2212,* 

4y 5*^oo,42i2: 



g •-'00,4212! 
^2212,* 



^00,4212! 



•-"00,4413 ~ •-'00,4211 ~ ^! 

^1112,* _ ^0011.* _ r, 
•-"00,5312 ~ •-'00,6211 ~ '^^ 



•'00,621 



(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 

(B18) 

(B19) 

(B20) 

(B21) 
(B22) 
(B23) 

(B24) 

(B25) 

(B26) 

(B27) 

(B28) 

(B29) 

(B30) 
(B31) 



with the two coupling constants CgQ qqqq left unrestricted. 
Apart from these 6 free and 36 dependent coupling con- 
stants, the vector-isoscalar channel of the CE requires 
that all the remaining sixth-order coupling constants are 
forced to be equal to zero. As before, we listed explicitely 
the vanishing coupling constants (see Eqs. (IB30I) - (|B31|) ). 
which were found to be non-vanishing in the scalar- 
isoscalar channel. 
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The general rule that we have specified at the end of 
the Appendix \^ can be applied now to explain the re- 
sults of this section, where at all the orders we found 
constraints that are nondiagonal in the spin space. Here, 
the reason is the possibility of having pairs of secondary 
densities that are coupled to rank v = 1. These pairs 
can appear at both scalar and vector coupling constants, 
which results in relating them to one another. 



At sixth order we have, 



Appendix C: Constraints for the vector-isovector 
channel (fourth and sixth orders) 



At fourth order we found the following constraints, 



0000, t 
00,4000 



= (^) 



^0011,0 
t '-'00,4011 

^0011,0 



r' 



2000,* 
00,2000 



1101,* 
00,3101 



r' 



2202,* 
00,2202 



3\/3 

.^0011,0 

v3 

{ /5„ooii,o 
3 V 3 00,40111 



1110,* 
00,3110 



a 



1111,* 

00,3111 



1112,* 
00,3112 



00,4011 



r' 



2011,* 
00,2011 



c, 



00,2211 



c, 



00,2212 



^^U01i,U 

-(^/3)*^V^C^^f, 
/:t^ooii,o 

V 'J*-^00,4011' 

^/^\t^r^ooii,o 

3*^00,4011' 
3^00,4011' 



5^0011,0 

3 '-^00,4011' 



c, 



2213,* 
00,2213 



(V3) 



C, 



00,2211 



3 
3 

1112.* 
00.3312 



7^0011,0 
3'-^00,4011' 

^0011,* 
'-'00,4211 



(CI) 

(C2) 

(C3) 

(C4) 
(C5) 

(C6) 

(C7) 
(C8) 
(C9) 

(CIO) 

(Gil) 

(C12) 
(C13) 



and only coupling constant C^q qJqq is unrestricted. 
Apart from these 2 free and 23 dependent coupling con- 
stants, the vector-isovector channel of the CE requires 
that all the remaining fourth-order coupling constants are 
forced to be equal to zero. In particular in the Eq. (|C13P 
we showed the vanishing coupling constants, which were 
found to be non- vanishing in the scalar-isoscalar channel. 



r 



■0000,* 
00,6000 



'^1 



2000,* 
00,4000 



•2202,* 
00,4202 



^1101,* 
•-"00,5101 



'^1 



3101,* 
00,3101 



'^1 



3303,* 
00,3303 



1110,* 
00,5110 



1111,* 

00,5111 



^1112,1 
'-'00,5112 



^3110,* 
•-"00,3110 



^3111,* 
'-'00,3111 



c; 



r" 



3112,* 
00,3112 

3312,* 
00,3312 



^3313,* 
'-'00,3313 



r" 



3314,* 
00,3314 



r' 



OOll,* 
00,6011 



■2011,* 
00,4011 



2211,* 
00,4211 



'^1 



2213,* 
00,4213 



^2212,* 
'-'00,4212 



,2211,* 
00,4011 

,2011,* 
00,4211 



-(V3)* 



C, 



1112.0 



2V15 

•,1112,0 



-(73)*-^^™^''"' 



2V15 

-,1112.0 



_(V3)*^^H|11£, 
( ^/3) * C'oo , 5if 2 ' 



1112,0 
00,51125 



(V3) 



^1112,0 
t '-'00,5112 



V5 



— >, V Y 1^^00,5112' 

y^lll2,0 

— '^'*^' 

= (^)*5y'^^00,5l'l2' 

— I fo\t'J_^lll2,0 

— \\ ii) g'-^00,5112' 

— I V o; g(-^00,5112' 

9^1112,0 

= (73)* 



= -(V3) 



3\/5 



00.5112 

y^lll2,0 

.(^)*7C^^ 

r,^1112,0 

_(y3)* ^^00,5112 



-(\/3)*2 



7 



C, 



15 

1112.0 



00,5112' 



-(V3)* ^^™ , 

^2213,* _ ^3112,* 
'-'00,4413 ~ '-'00,3312 



- r 



0011,* 
00,6211 



1112,* 
00,5312 

'0000. 
60,0000 



0, 
0, 



(C14) 

(C15) 

(C16) 

(C17) 

(C18) 

(C19) 

(C20) 

(C21) 
(C22) 

(C23) 

(C24) 
(C25) 
(C26) 

(C27) 
(C28) 
(C29) 
(C30) 
(C31) 
(C32) 

(C33) 

(C34) 
(C35) 



and only coupling constant Cqq qq'qq is unrestricted. 
Apart from 2 free and 39 dependent coupling constants, 
the vector-isovector channel of the CE requires that all 
the remaining sixth-order coupling constants are forced 



12 



to be equal to zero. In particular, in the Eqs. (jC34p - 
(|C35p we showed the vanishing coupling constants that 
were found to be non-vanishing in the scalar-isoscalar 
channel. 

The results presented in this section show simultane- 
ously both features we saw respectively in Appendices |X] 
and |B] At all orders, one can express all coupling con- 



stants through only one independent coupling constant, 
in such a way that the constraints are nondiagonal in 
both spin and isospin space. Again, this fact is due to 
the rank v = 1 and t = 1 in the pairs of densities in 
the final form of condition (|46)) . which allows the cou- 
pling constants at different spins and isospins to enter 
into the same constraints. 
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